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Abstract: - This paper studies two stochastic bulk arrival and bulk service C server queues (A) and (B) with k 
varying environments. The arrival and service times are exponential random variables and their parameters 
change when the environment changes. The system has infinite storing capacity and the arrival and service sizes 
are finite valued random variables. Matrix partitioning method is used to study the models. In Model (A) the 
maximum of the arrival sizes M in all the environments is greater than the maximum of the service sizes N in all 
the environments, (M > N), and the infinitesimal generator is partitioned as blocks of k times the maximum of 
the arrival sizes for analysis. In Model (B) the maximum of the arrival sizes M in all the environments is less 
than the maximum of the service sizes N in all the environments, (M < N), where the generator is partitioned 
using blocks of k times the maximum of the service sizes. Five different cases associated with C, M and N due 
to partitions are treated. They are namely, (Al) M >N > C, (A2) M >C >N (A3) C >M >N, which come up in 
Model (A); (Bl) N > C and (B2) C >N, which come up in Model (B) respectively. For the cases when C < M or 
N Matrix Geometric results are obtained and for the cases when C > both M and N Modified Matrix Geometric 
results are presented. The basic system generator is seen as a block circulant matrix in all the cases. The 
stationary queue length probabilities, its expected values, its variances and probabilities of empty queue levels 
are derived for the models using Matrix Methods. Numerical examples are presented for illustration. 
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I. INTRODUCTION 

In this paper two bulk arrival and bulk service multi server queues have been studied with random 
environment using matrix partitioning methods. Rama Ganesan, Ramshankar and Ramanarayanan in [1] and [2] 
have treated M/M/l bulk queue with random environment and PH/PH/1 bulk queue using Matrix Geometric 
Methods. Bini, Latouche and Meini [3] have studied numerical methods for Markov chains. Chakravarthy and 
Neuts [4] have discussed in depth a multi-server queue model. Gaver, lacobs and Latouche [5] have treated 
birth and death models with random environment. Latouche and Ramaswami [6] have studied Analytic 
methods. For matrix geometric methods and models one may refer Neuts [7]. The models considered in this 
paper are general compared to existing queue models since a multi server bulk arrival and bulk service queue 
with random environment has not been studied at any depth so far. Here random number of arrivals and random 
number of services are considered at a time whereas a fixed number of customers arrive or are served at any 
arrival or service epochs in many existing queue models in literature. Bulk service queue model, with service for 
fixed b customers when more than b customers are waiting, has been studied by Neuts and Nadarajan [8] for 
single server system. In the models considered here, depending on the environment the service and arrival sizes 
are random with distinct discrete distributions and the parameters of exponential distributions of arrival and 
service times vary. The number of servers increases with number of customers till it becomes C. The model with 
random arrival is similar to production systems. When a machine manufactures a fixed number of products in 
every production schedule, the defective items are always rejected in all productions; making any production lot 
is only of random size and not a fixed one always. Bulk service situations are seen often in software based 
industries where finished software projects waiting for marketing are sold in bulk sizes when there is economic 
boom and the business may be very dull when there is economic recession. In industrial productions, bulk types 
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are very common. Manufactured products arrive in bulk sizes and several bulk sizes of products are sold in 
markets. Recently M/M/l queue system with disaster has been studied by Noam Paz and Uri Yechali [9] but 
random arrival size or random service size with varying environments is not studied. Usually the partitions of 
the bulk arrival models have M/G/l upper-Heisenberg block matrix structure with zeros below the first sub 
diagonal. The decomposition of a Toeplitz sub matrix of the infinitesimal generator is required to find the 
stationary probability vector. Matrix geometric structures have not been noted so far as mentioned by William J. 
Stewart [10]. But in this paper the partitioning of the matrix is carried out in a way that the stationary probability 
vectors have a Matrix Geometric solution or a Modified Matrix Geometric solution for infinite capacity C server 
bulk arrival and bulk service queues with randomly varying environments. 

Two models (A) and (B) on M/M/C bulk queue systems under k varying environments with infinite 
storage space for customers are studied here using the block partitioning method. In the models considered here, 
the maximum arrival sizes and the maximum service sizes may be different for different environments. Model 
(A) presents the case when M, the maximum of all the maximum arrival sizes in all the environments is bigger 
than N, the maximum of all the maximum sale sizes in all the environments. In Model (B), its dual, N is bigger 
than M, is treated. In general in Queue models, the state space of the system has the first co-ordinate indicating 
the number of customers in the system but here the customers in the system are grouped and considered as 
members of M sized blocks of customers (when M >N) or N sized blocks of customers (when N > M) for 
finding the rate matrix. Using the maximum of the bulk arrival size or maximum of the bulk service size and 
grouping the customers as members of the blocks for the partitioning the matrix of the infinitesimal generator 
along with the environment state is a new approach in this area. In [1], Rama Ganesan, Ramshankar and 
Ramanarayanan for single server M/M/l bulk queue, have noticed two cases namely M >N and N >M but in this 
paper because a multi server system is of interest five cases are noticed. Model (A) gives three cases namely 
(Al) M > N > C, (A2) M > C > N and (A3) C > M > N and Model (B) gives two cases namely (Bl) N > C, and 
(B2) C > N. The case M=N with various C values can be treated using Model (A) or Model (B). For the cases 
when C < M or N, Matrix Geometric results are obtained and for the cases when C > both M and N, Modified 
Matrix Geometric results are presented. The matrices appearing as the basic system generators in these models 
due to block partitions are seen as block circulant matrices. The stationary probability of the number of 
customers waiting for service, the expected queue length, the variance and the probability of empty queue are 
derived for these models. Numerical cases are presented to illustrate their applications. The paper is organized in 
the following manner. In section II and section III the M/M/C bulk arrival and bulk service queues with 
randomly varying environment in which maximum arrival size M is greater than maximum service size N and 
the maximum arrival size is M less than the maximum service size N are studied respectively with their sub 
cases. In section IV numerical cases are presented. 

II. MODEL (A). MAXIMUM ARRIVAL SIZE M IS GREATER THAN MAXIMUM 

SERVICE SIZE N 

2.1 Assumptions for M > N. 

i) There are k environments. The environment changes as per changes in a continuous time Markov chain with 
infinitesimal generator Q 1 of order k with stationary probability vector ((). 

ii) Customers arrive in different bulk sizes for service. The time between consecutive bulk arrivals of customers 
has exponential distribution with parameter/!.; , in the environment i for 1 < i < k. At each bulk arrival in the 
environment i, Xi customers arrive with probability P (Xi= j) = pj for 1 < j < M t and Y^fli Pj=l for 1 < i < k. 

iii) Customers are served in batches of different bulk sizes. There are s servers to serve when s customers are 
present in the system for 1< s < C. When C or more than C customers are present in the system the number of 
servers to serve customers is C. In the environment i fori < i < k, the time between consecutive bulk services 
has exponential distribution with parameter sfx^ when s customers are in the system for 1< s < C and with 
parameter CfJ. t when C or more than C customers are present where jx t is the parameter of single server 
exponential service time distribution. At each service epoch in the environment i, ip t customers are served with 
probability given by P (tpi = j) = q) for 1< j < N t when more than N t customers are waiting for service 

where 'E^Ly q) =1- When n customers n < N t are in the system, then j customers are served with probability, qj 



International organization of Scientific Research 



14 I P a g e 



M/M/C Bulk Arrival And Bulk Service Queue With Randomly Varying Environment 



for 



1< j < n-1 and n customers are served with probability qj for 1 < i < k. 

iv) When the environment changes from i to j, the parameter of time between consecutive bulk arrivals and the 
service parameter change from (A t to (Aj,[ij ), the bulk arrival size Xi changes to Xj> the bulk service size 
ipi changes to ipj and the maximum arrival size M ; and the maximum service size N t change to Mj and Nj for 
l<i,j<k. 

v) The maximum of the maximum of arrival sizes M = max 1 <i<k M t is greater than the maximum of the 
maximum of service sizes N =max 1 <; <j< JV ; . 



2.2 Analysis 

There are three sub cases for this model namely (Al) M > N > C, (A2) M > C >N and (A3) C > M >N. Sub 
Cases (Al) and (A2) admit Matrix Geometric solutions and they are treated in sub section (2.2.1). Modified 
Matrix Geometric solution is presented for Sub Case (A3) which is studied in sub section (2.2.2). The state of 
the system of the continuous time Markov chain X (t) under consideration is presented as follows. 
X (t) = {(n, j, i): for 0 < j < M-l; 1 < i < k and n > 0} 

(1) The chain is in the state (n, j, i) when the number of customers in the system is n M + j, for 0 < j < M-l and 
0 < n < oo and the environment is i for 1 < i < k. When the number of customers in the system is r, then r is 
identified with (n, j) where r on division by M gives n as the quotient and j as the remainder. Let the survivor 
probability of the number of arrivals at an arrival epoch and the number of services at a service epoch 
in the environment i for 1 < i < k be P (jfj >j) = Pf =1-2^=1 Pn> and = 1 f° r 1 < j < M ; - 1 (2) 
and P M >j) = Q) =l-li =1 q l n , and Q l 0 = 1 for 1 <j <N t -1 (3) 
2.2.1 Sub Cases: (Al) M > N > C and (A2) M > C > N 

When M>N>CorM>C>N, the M/M/C bulk queue admits matrix geometric solution as follows. The chain 
X (t) describing them, has the infinitesimal generator Q A 21 of infinite order which can be presented in block 
partitioned form given below. 



<M,2.1 - 



B 1 A 0 0 

A 2 A 1 A 0 

0 A 2 A 1 

0 0 A 2 

0 0 0 



0 
0 



(4) 



In (4) the states of the matrices are listed lexicographically as 0, 1, 2, 3, ... . n, .... Here the vector n is of type 1 x 
k M and n = ((n, 0, l),(n, 0, 2)...(n, 0, k),(n, 1, l),(n, 1, 2)...(n, 1, k)...(n, M-l, l),(n, M-l, 2)... 
(n, M-l, k)) for n > 0. The matrices B 1 and A 1 have negative diagonal elements, they are of order Mk and their 
off diagonal elements are non- negative. The matrices A 0 ,andA 2 have nonnegative elements and are of order 
Mk and they are given below. Let the following be diagonal matrices of order k 

4=diag (hPj.hPj hvf) for 1 <j < M; U } = diag (^qj,n 2 qf l* k qf)for 1 < ; < N (5) 

Vj = diag (^Qj.^Qf H k Qf)forl <j < N; A=diag(A 1 ,A 2 A k ); U = diag (pi lt n 2 n k ) (6) 

let Qi = Qi - A - CU. (7) 
Here Q 1 is the infinitesimal generator of the Markov chain of the environment defined earlier 





0 


0 


0 


0 - 


Am-i 


A M ■ 


0 


0 


0 


Am-2 


Am-i " 


0 


0 


0 


Am-3 


Am-2 


0 


0 


0 


A3 


A, ■ 


•• A M 


0 


0 


A 2 


A3 ■ 


" A u _ x 


A M 


0 


i A, 


A 2 ■ 


" A M - 2 


Am-i 


A M - 



(8) 



0 • 


• 0 


CU N 


CI/jv-i ■ 


• cu 2 


cu 1 


0 ■ 


• 0 


0 


CU N ■ 


■ cu 3 


cu 2 


0 • 


• 0 


0 


0 


■ cu N 


cu N _, 


0 • 


• 0 


0 


0 


0 


CU N 


0 • 


• 0 


0 


0 


0 


0 


0 ■ 


• 0 


0 


0 


0 


0 



(9) 



International organization of Scientific Research 



15 I P a g e 



M/M/C Bulk Arrival And Bulk Service Queue With Randomly Varying Environment 



Qi 




A 2 ■ 


' A M _ N _ 


-2 


Am-n-i 


Am-n 


Am -2 


Am-i 


cu t 


Q'i 


A, ■ 


" A M _ N _ 


-3 


Am-N-2 


Am-n-i 


Am-3 


Am-2 


cu 2 


cu l 


Qi ■ 


" A M _ N _ 


-4 


Am-N-3 


Am-N-2 


Am -4 


Am-3 


CU N 


cu N _, 




■ Q'i 




A, 


A 2 


" A M _ N _ 2 


Am-n-i 


0 


CU N 


cu N _ t ■ 






Qi 


Ay 


■ A M _ N _ 3 


Am-n-2 


0 


0 


CU N ■ 


•• cu 2 




cu t 


Qi 


' A M _ N _ 4 


Am-N-3 


0 


0 


0 


■ cu N 






CU N _ 2 ■ 


■ Qi 


At 


0 


0 


0 


0 




CU N 






Qi 



(10) 



The matrix B l for Sub Case (Al) where N > C and Sub Case (A2) where C > N are given below in (1 1) and 
(12) respectively. For the case when C=N, the matrix^ may be written by placing C in place of N in the N-th 
block row in (12) and there after the multiplier of Uj is C. Let Q 1 j= 2i — A — jU for 0 < j < C and Q 1 c = Q[ 
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The basic generator 2,4 of the bulk queue, which is concerned with only the arrival and the service, is a matrix 
of order Mk given above in (13) where Q " A =A 0 + A 1 + A 2 . It is well known that a square matrix in which each 
row (after the first) has the elements of the previous row shifted cyclically one place right, is called a circulant 
matrix. It is very interesting to note that the matrix Q A is a block circulant matrix where each block matrix is 
rotated one block to the right relative to the preceding block partition. 
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The probability vector w of (13) gives, wQ A =0 and w.e = 1. (14) 
Since the probability vector of the environment generator Q 1 is <(>, the following are seen ^>Q 1 = 0 and (j) e =1. It 
can be seen in (13) that the first block-row of type k x Mk is W = (Q[ + A M ,A 1 ,A 2 , 
■ ■■A M _ N _ 2 , A M _ N _!, A M _ N + CU N , ... A M _ 2 + CU 2 , A M _ X + CU ± ) This gives as the sum of the blocks(2i + 

AM+ A1+ A2 + +AM-N-2+ AM-N-1 +AM-N+CUN+ +AM-2+CU2+ AM-1+CU1=Q1. Since 

<t>2i=0, this gives <|>(2i + A M ) + ifY.fJx N ~ x A t + * SUC^M-i + CU i) = 0 which implies (tp, <>...<(>,<(>). W = 0 = 
(()),())...((), (|) ) W where W is the transpose column vector of W. Since all blocks, in any block-row are seen 
somewhere in each and every column block due to block circulant structure, the above equation shows the left 

eigen vector of the matrix Q" A is {if, if ... if). Using (14) w = {m'm'm' -'m) (15) 

Neuts [7], gives the stability condition as, w A 0 e < w A 2 e where w is given by (15). Taking the sum of the 
same cross diagonally using the structure in (8) and (9) for the A 0 and A 2 matrices, it can be seen that 
wA 0 e =^(^ =1 nAje = ± f faE 0cDA 2 E iXi) hE(x k )) < wA 2 e = X -C $(Z% =1 nU n )e 

Taking the probability vector of the environment generator Q 1 as <|) = ,§ 2 , ... ,^ k ) , the inequality 

reduces to Sf=i <|> £ hE(xd < CH?=i § t HiE(&i )■ (16) 
This is the stability condition for the M/M/C bulk arrival, bulk service queue with random environment for 
Sub Case (Al) M > N > C and Sub Case (A2) M > C > N. When (16) is satisfied, the stationary distribution 
exists as proved in Neuts [7] . Let n (n, j, i), for 0 < j < M- 1 , 1 < i < k and 0 < n < oo be the stationary probability 
of the states in (1) and 7T n be the vector of type lxMk with, n n = (7i(n, 0, 1), 7i(n, 0, 2) ... 7i(n, 0, k), 7t(n, 1, 1), 
Ji(n, 1, 2),... 7t(n, 1, k)...7i(n, M-l, 1), 7i(n, M-l, 2)...7t(n, M-l, k) )for n > 0. The stationary probability 



vector 7T = (n 0 ,n 1 ,Ti 3 , ) satisfies 7r(?x i 2.i=0 and 7re=l. (17) 

From (17), it can be seen n 0 B 1 + 7T 1 v4 2 =0. (18) 

= 0, for H > 1. (19) 

Introducing the rate matrix R as the minimal non -negative solution of the non -linear matrix equation 

Ag+RA^^A^O, (20) 

it can be proved (Neuts [7]) that n n satisfies n n = n 0 R n for n > 1. (21) 

Using (18) and (21), n 0 satisfies n 0 [B r + RA 2 ] = 0 (22) 



Now 7T 0 can be calculated up to multiplicative constant by (22). From (17) and (21) n 0 (I — fi) _1 e =1. (23) 
Replacing the first column of the matrix multiplier of n 0 in equation (22) by the column vector multiplier of n 0 
in (23), a matrix which is invertible may be obtained. The first row of the inverse of that same matrix is n 0 and 
this gives along with (21) all the stationary probabilities. The matrix R given in (20) is computed using 
recurrence relation R(0) = 0; R(n + 1) = -A 0 Aj 1 -R 2 {n)A 2 A{ 1 , n > 0. (24) 
The iteration may be terminated to get a solution of R at an approximate level where \\R(n + 1) — R (n) || < e 
Note: 

The partition of the infinitesimal generator for the case M = C is similar and in that case C does not appear as a 
multiplier for the Uj matrices in the matrix B 1 (11) and (12) in the 0 block of (4) and C appears as a multiplier 
for all Uj matrices in the matrices of A 1 and A 2 from row block 1 on wards. From the arguments presented 
earlier it can be seen that the system admits Matrix Geometric solution for C = M also. 
2.2.2 Sub Case: (A3) C > M > N 
When C > M > N, the M/M/C bulk queue admits a modified matrix geometric solution as follows. The chain 
X (t) describing this Sub Case (A3), can be defined as in (1) presented for Sub Cases (Al) and (A2). It has the 
infinitesimal generator Q A 2 2 of infinite order which can be presented in block partitioned form given below. 
When C > M, let C = m* M + n* where m* is positive integer and n* is nonnegative integer with 0 < n* < M-l . 
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In (25) the states of the matrices are listed lexicographically as 0, 1, 2, 3, .... n, .... Here the vector n is of type 
1 xkMand n = ((n, 0, l),(n, 0, 2)...(n, 0, k),(n, 1, l),(n, l,2)...(n, 1, k)...(n, M-l, l),(n, M-l, 2)...(n, M-l, k)) 
for n > O.The matrices B \and A 1 have negative diagonal elements, they are of order Mk and their off diagonal 
elements are non- negative. The matrices A 0 andA 2 have nonnegative elements and are of order Mk and the 
matrices A^^-^and A 2 are same as defined earlier for Sub Cases (Al) and (A2) in equations (8), (9) and (10). 
Since C > M the number of servers in the system s equals the number of customers in the system L up to 
customer length C. When the number of customers L becomes more than C, (L > C), the number of servers in 
the system becomes constant C. When the number of customers L becomes less than C (L<C), the number of 
servers reduces and equals the number of customers. The matrix A 2 j for 1 < j < m*-l is given below 
0 ■•• 0 jMU N jMUfj^ ••• jMU 2 )MU 1 
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(26) 



The matrix A 2 mt is as follows given in (27) when C = m*M + n* and n* is such that 0 < n* < N-l. Here the 
multiplier of Uj in the row block increases by one till the multiplier becomes C = m*M + n* and there after the 
multiplier is C for Uj for all blocks. 

When N < n* < M-l, A 2 mt is same as in (26) for j = m*. 
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0 0 0 0 - 0 - 0 (Af-Wv - (Al-l)tfj CLw-j. 



Here^ ,= d - A - jU for 0 < j < C and Si c " Si-The matrix A! , for 1 <j < m*-l is as follows. 
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0 0 0 - 0 (JM + M-1)U K - (/Af+ Af-1)£/, C^x.K-.. 



The matrix mt is as follows when C = m*M + n* and n* is such that 0 < n* < N-l. The multiplier of Uj 
increases by one till it becomes C = m*M + n* and thereafter in all the blocks the multiplier of t/,- is C. 





4 


4 






4m 


4m 


(Aim. +l)i/ s 




4 




Air-pf-i 


Am-j 


4m 


(Aim. +?}£/, 






Am m < 




4m 


4m 



Cb',, 






• file 


4 


4i # i 


•Vv-.v-i 


0 


CI/ V 


CUm • 


CI/. 


Cu • 


4i w i 




0 


b 


0 


cos, 


CUm • 


■ <?U 


4 


0 




0 


0 


CU K 







When n* = N or n* > N then, in the matrix A 1 mt , there is slight change in the elements. When n* = N, in the 
N+l block row and thereafter C appears as multiplier of Uj , and when n* > N with n* = N + r for 1 < r < M-N-l, 
in the n*+l block row U N appears in the r + 1 column block. C appears as multiplier for it and as the multiplier 
of Uj thereafter in all row blocks respectively. The basic system generator for this Sub Case is same as (13) with 
probability vector as given in (15). The stability condition is as presented in (16). Once the stability condition is 
satisfied the stationary probability vector exists by Neuts [7]. As in the previous Sub Cases, 



7r( ?A,2.2=0 and 7re=l. (31) 
The following may be noted. TT n A 0 +n n+1 A 1 +n n+2 A 2 = 0, for n > m*, the rate mattix R is same as in previous 
Sub Cases with same iterative method for solving the same and n n satisfies n n = n mt R n ~ m * for n > m*. (32) 
The set of equations available from (31) are n 0 B' 1 +n 1 A 2 i= 0, (33) 

niA 0 +n i+1 A u+1 +n l+2 A 2ii+2 = 0, for 0 < i < m*-2 (34) 

and Tt m *-l A 0+ n m* A l,m*+ n m*+l A 2 = 0. (35) 

The equation 7re=l in (31) gives £™*o~ l7r ; e + 7r m*( i "R)~ le = 1 (36) 
Using Ti mt+1 = n mt R and equations (33), (34), (35) and (36) the following matrix equations can be seen. 
(n 0 ,n 1 ,n 3 , n mt ) Q' A22 = 0 (37) 

(n 0 ,n 1 ,n 3 , n mt ) _ R yi e ] =1 (38) 

The matrix Q' A22 is given by (39). 
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Q A,2.2~ 



B\ A 0 0 
A 2 .i A u A 0 
0 A 2i2 A u 



0 



0 



1 2,3 



0 
0 

^1,3 



0 
0 
0 

An 



0 
0 
0 
0 



(39) 



0 0 0 0 0 - A 2imt RA 2 +A lim 

Equations (37) and (38) may be used for finding (n 0 ,n 1 ,n 3 , n mt ). Replacing the first column of the first 

column- block in the matrix given by (39) by the column vector multiplier in (38) a matrix which is invertible 

can be obtained. The first row of the inverse matrix gives (n 0 ,n 1 ,n 3 , n mt ). 

This together with equation (32) gives all the probability vectors for this Sub Case. 
2.3. Performance Measures 

(1) The probability P(S = r), of the queue length S = r, can be seen as follows. Let m > 0 and n for 0 < n < M-l 
be non-negative integers such that r = mM+n. Then it is noted that 
P (S=r) =Hi=i n (. m > n > i )> where r = m M + n. 

(2) P (Queue length is 0) = P (S=0) = X?=i ^(0, 0, i). 

(3) The expected queue level E(S), can be calculated as follows. 
For Sub Cases (Al) and (A2) it may be seen as follows. Since n (n, j, i) = P [S = M n + j, and environment state 
= i], for n>0, and 0 < j < M-l and 1 < i < k, E(S) ^^oXf^ 1 If=i n(n,j, i) (Mn+j) 
=Zn=o n n- (Mn... Mn, Mn+1... Mn+1, Mn+2. . .Mn+2. . . Mn+M- 1 . . . Mn+M-1) where in the multiplier vector 
Mn appears k times, Mn+1 appears k times and so on and finally Mn+M-1 appears k times. 
So E(S) =MY l n= 0 nn n e +n 0 ( I - . Here Mk xl column vector §=(0, ...0,1. -,1,2, ...,2, ...,M - 
1,...,M—T where the numbers 0, 1, 2, 3,... and M-l appear k times in order. This gives 
E (S) = 7r 0 ( / - + Mn 0 (I - R)~ 2 Re (40) 
For Sub Case (A3), E(S) =E?= 0 £f B 5 l E*=i n(n,j, Q (Mn + ;) = M^ =0 nn n e + I,% =0 n n f = 
MSn =0 nn n e+ YiT=o 1 n A + n m*Q~R) Letting the generating function of probability vector O(s) — TiT=o n i s ' '> 
it can be seen, <D(s) = X™ q -1 n i si + n m* s m *(I-Rs) _1 and ^=0 nn n e = <D'(l)e = Y.?= *o _1 in t e+7r mt m*(I-R) _1 e 
+ 7r mt (I-R) -2 Re. This gives 

E(S) = M Er=o" 1 ^ i e + 7r mt m*(I-R)- 1 e + 7r mt (I-R)- 2 Re] +X™o~ 1 ^ + n mt (l-R)-^ (41) 

(4) Variance of queue level can be seen using Var (S) = E (S 2 ) - E(S) 2 . Let n be column vector n = 
[0,..0,l 2 ,...l 2 2 2 ,..2 2 ,... (M - l) 2 ,... (M - l) 2 ]' of type Mkxlwhere the squares of the numbers 0,1,2.. .(M- 
1) appear k times each in order. Then it can be seen that the second moment, for Sub Cases (Al) and (A2) 
E (S 2 ) = X =0 Y.fJo 1 EfLi n(n,j, i)[Mn+ j] 2 =M 2 [E» =1 n(n - l)n n e + E» =0 nn n e] + I» =0 + 
2M^ =0 nn n f. 

So, E(S 2 )=M 2 [n 0 (I - R)~ 3 2R 2 e + n 0 (I - Ry 2 Re]+n 0 (I - Ry 1 ^ + 2M n 0 Q - R)~ 2 R^ (42) 
Using (40) and (42) the variance can be written for Sub Cases (Al) and (A2). 
For the Sub Case (A3) the second moment can be seen as follows. E (S 2 ) = Xn=o Yif=o X?=i n ( n >ji 0 [M n + 
J]2= M2n=l<xmn- l7rne+n=(kon7rne+n=(ko7mrj +2Mn=(kon 7rn{ 

=M 2 [<D"(l)e +Zilo _1 to« e+7r mt m*(I-R)- 1 e + 7r mt (I-R)- 2 Re] + X^f 1 ^ + n mm (I-R)-\ 
+2ME™*o _1 in i f+7r mt m*(I-R)- 1 ^+7r mt (I-R)- 2 Ry. This gives E (S 2 ) = M 2 [ I]™*! -1 i(i - i)n t e + m*(m*- 
l)n mt (I - i?) _1 e +2m*?r mt (I-R)" 2 Re +27r mt (I-R)- 3 R 2 e +Y l fJ 0 ~ 1 ™i e+7r mt m*(I-R)- 1 e + 7r mt (I-R)- 2 Re] 
+ Sr^o" 1 n t r\ + ^.(I-R)" 1 !! +2M [ S^o" 1 in, f+7r mt m*(I-R)-^ + 7r mt (I-R)- 2 R Q. (43) 
Using (41) and (43) the variance can be written for Sub Case (A3). 



III. MODEL (B). MAXIMUM ARRIVAL SIZE M IS LESS 

THAN MAXIMUM SERVICE SIZE N 

In this Model (B) the dual case of Model (A), namely the case, M < N is treated. Here the partitioning matrices 
are of order Nk and the customers are considered as members of N blocks, M plays no role in the partition 
where as it played the major role in Model (A). Two Sub Cases namely (Bl) N > C and (B2) C > N come up in 
the Model (B). (When M =N and for various values of C greater than them, or less than them or equal to them, 
both models are applicable and one can use any one of them.) The assumption (v) of Model (A) alone is 
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modified without changing other assumptions stated earlier for the same. 
3.1 Assumption. 

v) The maximum of the maximums of arrival sizes in all the environments M = maxj <; </< M ; is less than the 
maximum of the maximums of service sizes in all the environments N =max 1 <; <j< JV ; where the maximum 
arrival and service sizes are M ; and jV ; in the environment i for 1 < i < k. 
3.2Analysis 

Since this model is dual, the analysis is similar to that of Model (A). The differences are noted below. The state 
space of the chain is as follows defined in a similar way presented for Model (A). 
X (t) = {(n, j, i): for 0 < j < N-l for 1 < i < k and 0 < n < 00} (44) 
The chain is in the state (n, j, i) when the number of customers in the queue is, n N + j, and the environment 
state is i for 0 < j < N-l, for 1 < i < k and 0 < n < 00. When the customers in the system is r then r is identified 
with (n, j) where r on division by N gives n as the quotient and j as the remainder. 
3.2.1 Sub Case: (Bl) N > C 

The infinitesimal generator Q B31 of the Sub Case (Bl) of Model (B) has the same block partitioned structure 
given in (4) for the Sub Cases (Al) and (A2) of Model (A) but the inner matrices are of different orders and 
elements. 



Qb,3.i- 



B\ 
A" 2 

0 

0 

0 



A\ 
0 
0 



0 

A\ 

A\ 
0 



0 
0 

A\ 



0 



0 



(45) 



In (45) the states of the matrices are listed lexicographically as 0, 1,2,3, .... n, ....Here the state vector is given 
as follows, n = ((n, 0, l),...(n, 0, k),(n, 1, l),...(n, l,k),(n, 2, l),...(n, 2, k),...(n, N-l, l),...(n, N-l, k)), for 0 < n 
< 00. The matrices, B'\, A" 0 ,A'\ and A" 2 are all of order Nk. The matrices B" 1 andA" 1 have negative 
diagonal elements and their off diagonal elements are non- negative. The matrices A" 0 andA" 2 have 
nonnegative elements. They are all given below. As in model (A), 
letting Aj ,, for 1 < j < M ,and Uj,Vj for 1 < j < N, A and U as diagonal matrices of order k given by (5) 
and (6) and letting Q[ = Q 1 — A — CU, the partitioning matrices are defined below 
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(46) 
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Here, Q' 1J = Q x - A -JU for 0 < j < C and Qj c = 2i- In (49) the case C > N has been presented. When C=N, 
Vj and in B'\ do not get C as multiplier in (49) and C appears as a multiplier of U j in A" 2 and A'\ in (47) 
and (48). The multiplier of matrices Uj and Vj concerning the services increase by one in each row block from 
third row block as the row number increases by one, up to the row C+l and it remains C in row blocks after that 
as given above. 
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The basic generator (50) which is concerned with only the arrival and service is Q B = A " 0 + A \ + A " 2 . This is 
also block circulant. Using similar arguments given for Model (A) it can be seen that its probability vector is 

W = ■"'«) an( ^ me stability condition remains the same. Following the arguments given for Sub Cases 

(Al) and (A2) of Model (A), one can find the stationary probability vector for Sub Case (Bl) of Model (B) also 
in matrix geometric form. All performance measures in section 2.3 including the expectation of customers 
waiting for service and its variance for Sub Cases (Al) and (A2) of Model (A) are valid for Sub Case (Bl) of 
Model (B) with M is replaced by N. It can also be seen that when N = C the system admits Matrix Geometric 
solution as in Model (A). 

3.2.2Sub Case: (B2) C > N 

The infinitesimal generator Qb$.2 of the Sub Case (B2) of Model (B) has the same block partitioned structure 
given in (25) for Sub Case (A3) of Model (A) but the inner matrices are of different orders and elements. When 
C > N > M, the M/M/C bulk queue admits a modified matrix geometric solution as follows. The chain X (t) 
describing this Sub Case (B2), can be defined as in the Sub Case (Bl). It has the infinitesimal generator Q B 3 2 of 
infinite order which can be presented in block partitioned form given below. When C > N, let C = m* N + n* 
where m* is positive integer and n* is nonnegative integer with 0 < n* <N-1. 
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In (51) the states of the matrices are listed lexicographically as 0, 1,2,3, ... .n, .... Here the vector n is of type 
1 x kN and n = ((n, 0, l),(n, 0, 2)...(n, 0, k),(n, 1, l),(n, 1, 2)...(n, 1, k)...(n, N-l, l),(n, N-l, 2)...(n, N-l, k) 
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for n > 0. The matrices B"\, A'\i for 1 < j < m* and A'\ have negative diagonal elements, they are of order 
Nk and their off diagonal elements are non- negative. The matrices A" 0 ,A" 2 j and A" 2 for 1 < j < m* have 
nonnegative elements and are of order Nk and the matrices A" 0 ,A" 1 and A" 2 are same as defined earlier for Sub 
Case (Bl) in equations (46), (47) and (48). Since C > N the number of servers in the system s equals the number 
of customers in the system L up to customer length becomes C= m* N + n*. When the number of customers 
becomes more than C, (L > C), the number of servers in the system becomes constant C. When the number of 
customers becomes less than C, the number of servers again falls and equals the number of customers. As in 
model (A), letting A, for 1 <j < M ,and Uj,V for 1 < j < N, A and U as diagonal matrices of order k 
given by (5) and (6) and letting Q[ = Q 1 — A — CU, the partitioning matrices are defined as follows. The matrix 
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(52) 



The matrix A 2 mt is as follows : 
(Nm*)U N (Nm*)U N _ 1 
0 (Nm * +1)U N 



iven in (53) when C = m*N + n* is such that 0 < n* < N-l. A " 2 _, 
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Considering again, ,= — A — jU for 0 < j < C and Q\ : = Si 33 m Sub Cases [Al) and(A2), 
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The matrix A'\ , for 1 < j < m*-l is as follows. 



QX+ IV - M - l}£/v-x-i C/' v + W - M - Wv-x-s 



G-.v-.v-2;::V- : 



0A-+ AT- l)£/ w -x Wv-,-x 



The matrix A" l mt is in (56) when C = m*N + n* and 0 < n* < N-l. From row block n*+l, the multiplier of Uj is 



C. The matrix A" lm , 
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The basic generator for this model is also same as (50) which is concerned with only the arrival and the service. 
Q B = A" 0 + A'\ + A" 2 . This is also block circulant. Using similar arguments given for Model (A) it can be 

seen that its probability vector is w = (~>Jj>~>—>Jj) an( ^ me stability condition remains the same. Following 
the arguments given for Sub Case (A3) in section 2.2.2 of Model (A), one can find the stationary probability 
vector for Sub Case (B2) of Model (B) also in modified matrix geometric form. All the performance measures 
given in section 2.3 including the expectation of customers waiting for service and its variance for Sub Case 
(A3) are valid for Sub Case (B2) of Model (B) except M is replaced by N. 



IV. NUMERICAL ILLUSTRATION 

For the models the varying environment is considered to be governed by the Matrix Q 1 = 
3 

2 . The arrival time and service time parameters of exponential distributions are respectively 



fixed in the three environments El, E2 and E3 as A = (10, 12, lA)and fx = (4,5,6). Twelve examples are 

studied with various values for C, M and N. The maximum arrival size M and the maximum service size N in all 
environments are (i) M=N=4 in four examples, (ii) M=4, N=3 in four examples studied for Model (A) and (iii) 
M=3, N=4 in four examples studied for Model (B). In all these sets of different values of M and N, mentioned in 
(i), (ii) and (iii) the number of servers C are varied as C = 3, 4, 6 and 7. When M=N=4 the probabilities of bulk 
arrival and bulk service sizes in the environments are as follows given in Table 1. 



Table 1: Probabilities of / Arrival and \\> Service Sizes with Maximums = 4 in Three Environments. 
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For the cases with M=3, without changing the arrival size probabilities given above in Table 1 for the 
environments E2 and E3, the arrival size probabilities in El are changed as follows P( Arrival size =1 ) =.5 
P( Arrival size =2) =.4, P( Arrival size=3)= .1 and P( Arrival size =4) =0. For the cases with N=3, without 
changing the service size probabilities given above in Table 1 for the environments E2 and E3, the service size 
probabilities in El are changed as follows P(Service size =1 ) =.5 P(Service size =2) =.3, P(Service size =3)= .2 
and P(Service size = 4) =0. Here same numbers of 30 iterations are performed to find the rate matrix R in all the 
twelve cases. When C =3 and C =4 matrix geometric results are seen and the results obtained are presented in 
Table2. When C= 6 and C=7 modified matrix geometric results are seen and the results obtained are presented 
in Table 3. Probabilities of various sizes of queue lengths S = 0, 1, 2, 3 and various blocks 0, 1, 2, 3 are 
obtained. Further P(S> 15), E(S) and VAR (S) are also derived for all the twelve cases. Norms, arrival rate and 
service rate values are presented. Figure 1 and figure 2 present probabilities of various queue sizes and block 
sizes for the twelve examples. Effects of variation of rates on expected queue length and on probabilities of 
queue lengths are exhibited in tables 2 and 3 respectively. The decrease in arrival rates (so also increase in 
service rates) makes the convergence of R matrix faster which can be seen in the decrease of norm values. 



Table2: Results Obtained For Matrix Geometric Models. 





C=3,M=4=N 


C=3,M=4,N=3 


C=3,M=3,N=4 


C=4=M=N 


c=4,M=4,N=3 


C=4,M=3,N=4 


P(S=0) 


0.079684166 


0.072645296 


0.086694481 


0.134207116 


0.127898311 


0.138717179 


P(S=1) 


0.094816229 


0.089299727 


0.103026395 


0.160159184 


0.158139375 


0.165224098 


P(S=2) 


0.099972009 


0.094739563 


0.10844421 


0.168932101 


0.167596052 


0.173956131 


P(S=3) 


0.080822795 


0.077354801 


0.088540336 


0.1362632 


0.136279652 


0.141667691 


7t0e 


0.3552952 


0.334039388 


0.386705422 


0.599561601 


0.589913389 


0.619565098 


nle 


0.246062643 


0.239774491 


0.258608447 


0.269420604 


0.27246095 


0.268090953 


n2e 


0.151711966 


0.152962345 


0.149421908 


0.087788452 


0.091032247 


0.079174934 
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7r3e 


0.093923215 


0.098003396 


0.086467581 


0.028918419 


0.030765921 


0.023368742 


P(S>15) 


0.153006976 


0.175220381 


0.118796643 


0.014310924 


0.015827493 


0.009800274 


Norm 


8.607020E-05 


0.000121229 


5.72322E-05 


1.70312E-07 


2.66577E-07 


6.89660E-08 
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4.322222222 


4.322222222 


4.22962963 


4.322222222 


4.322222222 


4.22962963 


Service 
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5.35 


5.238888889 


5.35 


7.133333333 


6.985185185 


7.133333333 


E(S) 


8.177499726 


8.832269818 


7.221716529 


3.754933893 


3.853312418 


3.52779164 


Var(S) 


69.84639946 


81.04471283 


53.74109494 


13.38646744 


13.99187868 


11.37540035 



Table3: Results Obtained For Modified Matrix Geometric Models 
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Figure 1: Probabilities of Queue lengths 



Figure 2: Probabilities of Block Queue Sizes 



V. CONCLUSION 

Two M/M/C bulk arrival and bulk service queues and their sub cases with randomly varying 
environments have been treated. The environment changes the arrival rates, the service rates, and the 
probabilities of sizes of bulk arrivals and bulk services. Matrix geometric and modified matrix geometric results 
have been obtained by suitably partitioning the infinitesimal generator by grouping of customers and 
environments together respectively when the number of servers is not greater than or greater than the maximum 
of the maximum arrival and maximum service sizes. The basic system generators of the queues are block 
circulant matrices which are explicitly presenting the stability condition in standard form. Numerical results for 
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various bulk queue models are presented and discussed. Effects of variation of rates on expected queue length 
and on probabilities of queue lengths are exhibited. The decrease in arrival rates (so also increase in service 
rates) makes the convergence of R matrix faster which can be seen in the decrease of norm values. The 
variances also decrease. Bulk PH/PH/C queue with randomly varying environments causing changes in sizes of 
the PH phases may produce further results if studied since PH/PH/C queue is a most general form almost 
equivalent to G/G/C queue. 
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